Relaxation and dephasing in a flux qubit 
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We report detailed measurements of the relaxation and dephasing time in a flux-qubit measured by 
a switching DC SQUID. We studied their dependence on the two important circuit bias parameters : 
the externally applied magnetic flux and the bias current through the SQUID in two samples. We 
demonstrate two complementary strategies to protect the qubit from these decoherence sources. 
One consists in biasing the qubit so that its resonance frequency is stationary with respect to the 
control parameters {optimal point) ; the second consists in decoupling the qubit from current noise 
by chosing a proper bias current through the SQUID. At the decoupled optimal point, we measured 
long spin-echo decay times of up to 
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PACS numbers: 

A long-standing problem for the use of superconduct- 
ing circuits as quantum bits (qubits) in a quantum com- 
puter Q, 1^ H, |j{ is their relatively short dephasing time 
compared to the requirements of many-qubit quantum 
computation. Dephasing is due to the coupling of the 
qubit 's degrees of freedom with the many fluctuating un- 
controlled ones commonly denoted as the environment 
P, 13 . From the perspective of quantum information, it 
is crucial to quantitatively identify the various dephasing 
sources and to find strategies to overcome these, either by 
reducing the amount of fluctuations or by protecting the 
qubit against it. An inxportant step in this direction has 
been accomplished in Q. The authors showed that de- 
phasing can be significantly reduced by biasing the qubit 
at an optimal point where its resonance frequency is sta- 
tionary with respect to its control parameters - in that 
case, gate voltage and magnetic flux. 

In this letter we report detailed measurements of the 
relaxation and dephasing times in a flux-qubit as a func- 
tion of its bias parameters for two different samples. 
Our measurements allow us to identify certain dephas- 
ing mechanisms and quantify their effect on the qubit. 
We find that energy relaxation is dominated by sponta- 
neous emission towards the measuring circuit impedance. 
Dephasing is mainly caused by noise in the external mag- 
netic flux biasing the qubit, thermal fluctuations of our 
measuring circuit, and low- frequency noise originating 
from microscopic degrees of freedom, probably causing 
critical current noise in the qubit junctions. We moreover 
demonstrate strategies to efficiently fight each of these 
noise sources. 

Our flux-qubit consists of a micron-size superconduct- 
ing loop intersected with three Josephson junctions @. 
When the total phase across the three junctions jq 
is close to TT, the loop has two low-energy eigenstates 
(ground state |0) and excited state |1)) well separated 
from the higher-energy ones, which can thus be used 
as a qubit |j, 01- The flux-qubit is characterized by 
two parameters : the minimum energy separation A be- 
tween |0) and and the persistent current Ip. Around 



7q = TT, the energy separation between these two levels 
depends on jq and can be written as Ei — Eq = hfg = 
e^, where e= {Ip / e) {jq — tt) / {2Tr) . The qubit is 



inductively coupled to a SQUID detector (with a coupling 
constant Af), which is biased at a current lb. The phase 
drop 7q has two origins : the magnetic flux threading the 
qubit loop (f>2,, and the currents in the SQUID loop which 
depend on Thus, we can write e = 'i]{^x) + ^{h)- 

The coupling of e to fluctuating sources leads to de- 
coherence. Noise in the magnetic flux or in the bias 
current induces fluctuations of the qubit frequency /q 
and thus dephasing. A flrst strategy to protect the qubit 
from decoherence consists in biasing it at e = so that 
d/g/rfe ~ 0. This is the optimal point strategy, which 
was first invented and demonstrated in 0| . An additional 
possibility is to decouple the external noise from the vari- 
able e, by canceling the sensitivity coefficients d-q/d^x 
and dX/dlb. The flux noise can not be decoupled since 
drj/d^x = 2Ip/h is constant. As we will show below, 
the bias current noise can be decoupled by biasing the 
SQUID at a current I^ such that dX/dIb{Ib) — 0, which 
is the decoupling condition. At the decoupled optimal 
point, (e = and If, = I^) we expect that the qubit 
quantum coherence is best preserved, since the qubit is 
sensitive to flux noise to second order, and to bias current 
noise to fourth order. We also note that a strong depen- 
dence of the dephasing time on the bias current would be 
clear experimental evidence that current noise, and not 
flux noise, is the factor limiting the quantum coherence. 

In the two samples, shown in figures ^ and the 
qubit loop is merged with its measuring SQUID. The de- 
pendence of e on the bias current If, arises from the way 
this bias current redistributes in the SQUID and even- 
tually generates a phase shift across the qubit junctions 
via the superconducting line shared by the qubit and the 
SQUID. The detailed configuration of the shared line is 
related to the specific fabrication process. We use 2-angle 
shadow evaporation so that the lines consist effectively of 
2 layers. This induces a large asymmetry in the coupling 
^ if the qubit loop contains an odd number of junc- 
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FIG. 1: (a) Atomic Force Micrograph of sample A. The flux 
qubit is the small loop containing three Josephson junctions 
in a row ; the SQUID is constituted by the outer loop contain- 
ing the two large junctions. The bar indicates a Ifim length, 
(b) Scanning Electron Micrograph of sample B. Note that 
the qubit loop contains a fourth junction, 3 times larger than 
the other ones, (c) Electrical model of the measuring circuit. 
The SQUID, represented by its Josephson inductance Lj, is 
shunted by an on-chip capacitor Csh through superconduct- 
ing lines of inductance L (all on-chip). It is current-biased 
by a waveform generator delivering a voltage E across an 
impedance Zin ; the voltage across the SQUID is connected 
to the input of a room-temperature preamplifler through an 
impedance Zout- Zin and Zout include on-chip gold resistors. 



tions. In this article, we compare the results obtained 
for a three- (sample A) and a four- (sample B) junction 
qubit (see figure^ and b), and demonstrate that such 
asymmetry can be removed by using an even number 
of junctions in the qubit loop A model for the qubit 
electromagnetic environment in both samples is shown in 
figure Ic. The SQUID is modeled by its Josephson induc- 
tance Lj shunted by a capacitor Csh via superconducting 
lines of inductance L. It is connected to the output volt- 
age of our waveform generator E via an impedance Zin, 
and to the input of a room-temperature amplifier through 
an impedance Zout- It thus forms a harmonic oscillator, 
the plasma mode, of frequency ujp = (a/ {L + Lj)Csh)~^ 
and quality factor Q = u)j^shRe-{Z){ujp), to which the 
qubit is strongly coupled 9] (we note Z — Zin/ / Zout)- 
Here is a list of the parameters for our two samples : 
for sample A, Ip = 270nA, A = 5.85GHz, M = 20pH, 
Lj = 80pH, L = llOpH, C = 12pF, Z{0) = lAkn ; 
for sample B, Ip = 240nA, A = 5.5GHz, M = 6.5piJ, 
Lj = 380pff, L = 80pff, Lj = 350pH, C = 5.5pF, 
Z{Q) = 9/cri. We note that the main difference between 
sample A and B, apart from the number of junctions in 
the qubit loop, is the value of M and of the low-frequency 
impedance Z{Q). 

We measured our sample parameters and A(/f,) by 
studying the dependence of the qubit Larmor frequency 
on both the external flux and the bias current 1^. We 
performed spectroscopy by applying a 500ns microwave 
pulse of variable frequency, and measuring the SQUID 
switching probability with a short subsequent DC cur- 
rent pulse |j] for different values of We added a l/xs 
plateau at the value hpi in order to adjust the bias cur- 
rent through the SQUID during the application of the 
microwave pulse. The complete pulse sequence is de- 
picted in figure 12^" In figure the measured qubit res- 
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FIG. 2: (a) Principle of the spectroscopy experiments : a 
bias current pulse of amplitude hbpl (lower than the SQUID 
critical current) is applied to the sample while a microwave 
(MW) pulse probes the qubit resonance frequency. The qubit 
state is finally measured by a short bias current pulse as 
discussed in j4|. (b) Typical spectroscopy curves for three 
values of I^pi measured with sample A (from left to right, 
hpi = — 2.25, 0, 2/^j4). The solid curves are fits to the data. 
(c,e) Curves \{Ib) deduced from the spectroscopy curves as 
explained from the text for sample A (c) and sample B (e). 
The decoupling condition is satisified at = 2.9 ± OA/iA for 
sample A (black arrow in the figure) and = 180 ± 20nA for 
sample B. (d,f) Qubit line at the decoupled optimal point for 
sample A (d) and B (f). 



onance frequency for sample A is shown as a function of 
the external flux for three different values of Ibpi ■ We 
observe that for each value of the bias current, a speciflc 

value of external flux ^i^'' (Itpi) realizes the optimal point 
condition. 

We fitted all the curves with the formula fg = 
v/A2 + [A(/;,) + 2/p($^ - $o/2)//i]2 for different values 
of /(,. The obtained curves A(Jf,) are shown in figure 
13; and 1213 for both samples. The decoupling occurs at 
4* = 2.9±0AfiA for sample A and at i; = 180±20nA for 
sample B. Note that although the SQUID critical current 
is similar in both samples, the decoupling current is much 
closer to in sample B due to the presence of the fourth 
junction wich restores the symmetry of the coupHng 0. 
We biased our qubit at the decoupled optimal point by 
setting lb — 4* and = ^l^Hlb)- The qubit line shape 
under these conditions is shown in figure for sam- 
ple A and|2t for sample B. For sample A, we could fit it 
with a Lorentzian of width w = 3A±0.5MHz (FWHM). 
This width yields a dephasing time T2 — l/nw ~ lOOns 
consistent with the Ramsey fringe measurements as dis- 
cussed below. For sample B, the line was split, due to 
the action of a strongly coupled two-level fiuctuator. We 
fitted it by the sum of two Lorentzians of widths 7 and 
6AIHz. We note that in addition to the fluctuator re- 
sponsible for the splitting of the line, the value of the 
qubit frequency at the optimal point A exhibited occa- 
sional jumps of around lOOMHz. Also the width of the 
line changed significantly in time. This indicates that de- 
phasing was probably dominated by some low-frequency 
noise due to one or more strongly coupled microscopic 
fluctuators, likely generating critical current noise. We 
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FIG. 3: (a) Ramsey fringe signal measured with sample A. 
From top to bottom : - temporal sequence of microwave pulses 
corresponding to a Ramsey experiment. - Ramsey signal at 
the optimal decoupled point. - Qubit resonance frequency as 
a function of e. - Dephasing time T2 as a function of e around 
the optimal point (full squares, the lines are a guide to the 
eye), and fit to the data (dotted curve) assuming that dephas- 
ing was caused by 1// flux noise = 3 • 10~^^//[$o/^2]. 
(b) Spin-echo signal measured with sample B. From top to 
bottom : Temporal sequence of microwave pulses correspond- 
ing to a spin-echo experiment. - Spin-echo signal at the opti- 
mal decoupled point. - Qubit resonant frequency as a function 
of e. - (full squares) : Spin-echo time Techo as a function of e 
(the lines are a guide to the eye) . (dashed curve) : calculated 
dephasing from thermal fluctuations of the photon number in 
the SQUID plasma mode. 



stress that we had no evidence for such instabihties with 
sample A. 

We first studied the dependence of the dephasing time 
as a function of e while keeping 7^ — I^- For sample A, we 
measured Ramsey fringes 3, 4] by applying a sequence 
of microwave pulses as schematized in figure for each 
value of <J>a;. The Ramsey fringes measured at the decou- 
pled optimal point are shown. They decay exponentially 
with a time constant T2. Figure (bottom) shows the 
dependence of T2 with $j; . The dephasing time exhibits 
a sharp maximum T2 = 120ns at the optimal point e = 
for which d/g/de = as expected. 

To account for the rapid degradation of the dephas- 
ing time when e ^ 0, we first evaluated the effect of 
the thermal fluctuations in the measuring circuit on the 
qubit coherence time. Thermal fluctuations of the pho- 
ton number in the plasma mode cause fluctuations of 



the qubit resonance frequency and thus dephasing. It 
has been shown that treating the thermal fluctua- 
tions of the plasma mode as a weak classical perturba- 
tion leads to a strong underestimate because of the ne- 
glect of quantum correlations between the qubit and the 
oscillator. Instead, we numerically integrated the mas- 
ter equation for the joint density matrix of the "qubit- 
plasma mode" system [l5| . The qubit density matrix is 
obtained at the end of the calculation by tracing over 
the plasma mode degrees of freedom. The evolution of 
its off-diagonal element yields the dephasing time. For 
these calculations, we assumed a quality factor Q — 100 
for sample A and Q = 150 for sample B, and an effective 
temperature T = 70mK in agreement with additional 
measurements '9^ ; all the other parameters of the model 
are directly obtained from experimental data. For sample 
A, we found that thermal fluctuations have a negligible 
effect when 1^ = 1^; we thus believe that flux-noise is 
responsible for rapid degradation of the dephasing time 
when e ^ 0. Assuming that the flux-noise power spec- 
trum has a frequency dependence given by S<s,^ = A/\f\ 
which is consistent with noise measurements found in the 
literature use the data fror n fig ure to eval- 

uate A. With calculations similar to [lj|, we find that 
A = 3± 1.5 • 10~^^<i>g gives a good agreement (dotted hue 
in figure Ell). Such a level of noise is comparable to the 
lowest values reported in SQUID measurements 0, ■ 

As could be expected from the lineshape shown in fig- 
ure|2f, the Ramsey fringe signal measured with sample B 
had a non-exponential damping so that it was impossi- 
ble to measure T2. To circumvent the low- frequency noise 
mentioned above, we used a spin-echo type sequence of 
microwave pulses shown in figure ^jp (top) , as demon- 
strated in the case of low-frequency charge noise 0, ^3 ■ 
The results are shown in figure |3b at the decoupled op- 
timal point, by a set of curves corresponding to different 
delays between the two tt/2 pulses. We fitted each curve 
by a gaussian multiplied by a sine curve. Then we fit- 
ted the decay of the echo amplitude as a function of the 
delay between the two 7r/2 pulses with an exponential of 
time constant T^cho- At the decoupled optimal point, we 
measured T^cho = 3.9 ± O.l^s. We stress that these re- 
sults represent a significant improvement over previously 
reported coherence times in superconducting qubits. We 
studied the dependence of Techo as a function of e (figure 
OId bottom, full squares) for If, = I^. Again we found a 
sharp maximum at e = 0. For sample B, thermal fluctu- 
ations in the plasma mode account qualitatively for the 
experimental data (dashed curve in figure bottom). 

We finally studied the bias current lb dependence of 
the dephasing and echo times T2 and T^cho together with 
the energy relaxation time Ti at e = 0. The results are 
shown in figure^Ji and c for samples A and B respectively. 
All these curves exhibit a clear maximum at /& = I^. 
This indicates that at If, ^ both relaxation and de- 
phasing are limited by coupling to the measuring circuit. 
In particular, the fact that in sample B Ti is strongly 
reduced away from is a clear indication that energy 
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FIG. 4: (a) Relaxation (Ti, black squares) and dephasing {T2, 
empty circles) times at the optimal point (e = 0) as a func- 
tion of the bias current for sample A. (b) Pure dephasing 
time (full squares) as a function of /;,. The dashed line 
is the result of a simulation taking into account the thermal 
fluctuations of the plasma mode, (c) Relaxation (Ti, black 
squares) and echo {T2, empty circles) times at the optimal 
point (e = 0) as a function of the bias current /(, for sample 
B. (d) Pure dephasing component of the echo time T^'^'^° (full 
squares) and calculated effect of the thermal fluctuations in 
the plasma mode (dashed line). 

relaxation occurs by spontaneous emission towards the 
circuit impedance seen by the qubit. A weaker depen- 
dence is observed for sample A, which could indicate 
that additional environmental modes at the qubit fre- 
quency are involved. The dephasing time T2 measured 
in sample A is strongly dependent on /;,. This indicates 
that dephasing at the optimal point is limited by noise in 
the bias current. For both samples, the dephasing time 



measured at the optimal decoupled point is similar to or 
larger (sample A) than the relaxation time, so that de- 
phasing was partly limited by relaxation. To quantify the 
pure dephasing contribution, we calculated defined as 
T^^ = (Ta)"^ - (2Ti)-i for sample A (fuU square curve 
in figure ^ja) and calculated similarly T^^^° for sample 
B (full square curve in figure^). Our calculations tak- 
ing into account the thermal fluctuations in the plasma 
mode are shown as the dashed curve in figures ^ and 
d. They are in qualitative agreement with the data, al- 
though systematically overestimating the dephasing time 
by a factor typically 5 compared to the measurements. 

In conclusion, we presented detailed measurements of 
the relaxation and dephasing times as a function of bias 
parameters for two flux-qubit samples. We showed that 
the optimal point concept already demonstrated for the 
quantronium circuit j^] is also valid for the flux-qubit de- 
sign. Making use of the SQUID geometry of our detector, 
we could moreover decouple the qubit from current fluc- 
tuations by biasing the SQUID at a speciflc current I^. 
We showed that adding a fourth junction to the qubit 
loop enhances the symmetry of the coupling, thus lower- 
ing the value of . We showed that low- frequency noise 
limits the dephasing time, but that spin-echo techniques 
provide a powerful tool to fight it. We observed remark- 
ably long decay times of the echo signal of 4/xs, limited 
by relaxation. We provided quantitative evidence that at 
the optimal point dephasing is induced by the thermal 
fluctuations of the photon number in the plasma mode 
of our SQUID detector. These results indicate that long 
coherence times can be achieved with flux qubits. 
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